Ferrimagnetism of the magnetoelectric compound Cu 2 OSe0 3 probed by 77 Se NMR 
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We present a thorough 77 Se NMR study of a single crystal of the magnetoelectric compound Cu20SeC>3. 
The temperature dependence of the local electronic moments extracted from the NMR data is fully consistent 
with a magnetic phase transition from the high-T paramagnetic phase to a low-T ferrimagnetic state with 3/4 of 
the Cu 2+ ions aligned parallel and 1/4 aligned antiparallel to the applied field of 14.09 T. The transition to this 
3up-ldown magnetic state is not accompanied by any splitting of the NMR lines or any abrupt modification in 
their broadening, hence there is no observable reduction of the crystalline symmetry from its high-T cubic />2i3 
space group. These results are in agreement with high resolution x-ray diffraction and magnetization data on 
powder samples reported previously by Bos et al. [Phys. Rev. B, 78, 094416 (2008)]. We also develop a mean 
field theory description of the problem based on a microscopic spin Hamiltonian with one antiferromagnetic 
(Jafm — 68 K) and one ferromagnetic (Jf m ~ —50 K) nearest-neighbor exchange interaction. 

PACS numbers: 76.60.Jx, 75.50.Gg, 75.25.-j, 77.84.Bw 



I. INTRODUCTION 

Multiferroic and magnetoelectric materials are currently at 
the center of intense research activity*^ In magnetoelectric 
compounds the application of an electric field can induce a fi- 
nite magnetization and similarly an electric polarization can 
be induced by an applied magnetic field.^ Combining elec- 
tronic and magnetic properties is an exciting playground for 
fundamental research which may also lead to new multifunc- 
tional materials with potential technological applications^iS 
Since both time-reversal and spatial-inversion symmetries 
must be broken in ferroic materials^ magnetoelectric (ME) 
effects are allowed only in 58 out of the 122 magnetic point 
groups^i^ Along with the weakness of the ME effect, this re- 
sults in a limited number of compounds displaying ME prop- 
erties, such as Cr 2 0;A Gd 2 CuO^, and BaMnF^. Even 
though ME materials have been studied extensively in the past 
decades the recent discovery of the multiferroic compounds 
TbMnCV 14 and TbMn 2 Og™, where ferroelectricity (FE) is 
driven directly by the spin order, has led to a revival of interest 
in ME systems. The magnetic state in these compounds is ei- 
ther a spiral (as in TbMnOg^ Ni 3 V 2 8 ^ and MnWO^iiH), 
or a collinear configuration (as in Ca3(CoMn) 2 C<6r2 and 
FeTe 2 5 Bs22). 

Considerable interest has also been drawn to the investi- 
gation of ME effects appearing in non-polar systems such as 
SeCuOs^ TeCuOa^ and CuaOSeO^. The latter, as re- 
ported recently by Bos et al.r^, undergoes a ferrimagnetic 
phase transition at T c ~ 60 K which is accompanied by a 
significant magnetocapacitance signal and an anomaly of the 
dielectric constant. However, high resolution powder x-ray 
diffraction (XRD) data show that the lattice remains metri- 
cally cubic down to 10 K, and this excludes a ME coupling 
mechanism that involves a spontaneous lattice strain. This 
result is further supported by recent infrared,™ and Raman 24 
studies. In this respect, Cu 2 OSe03 appears to be a unique 
example of a metrically cubic material that allows for piezo- 
electric as well as linear magnetoelectric and piezomagnetic 




FIG. 1: (Color online) Crystal structure of Cu20Se03. Copper atoms 
(Cui in blue and Cu2 in green) form a network of distorted corner 
sharing tetrahedra. Only the ions inside the unit cell are shown. Oxy- 
gen atoms are shown in yellow while Sei and Sen ions are shown in 
magenta and red respectively. 



coupling. 

Here we present an extensive 77 Se Nuclear Magnetic Res- 
onance (NMR) study which highlights the local magnetic and 
structural properties of Cu 2 OSeC<3. Our NMR measurements 
are performed in a single crystal of Cu 2 OSe03 as a function of 
temperature and by varying the direction of the applied mag- 
netic field with respect to the crystalline axes. A detailed anal- 
ysis of the 77 Se spectra based on space group symmetry con- 
siderations leads to the following main conclusions: The tem- 
perature dependence of the local electronic moments extracted 
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from the NMR data is fully consistent with a phase transition 
between the high-T paramagnetic phase and a low-T ferrimag- 
netic configuration with 3/4 of the Cu 2+ moments (type CU2) 
aligned parallel and 1/4 (type Cui) aligned antiparallel to the 
applied field. This 3up-ldown state has 1/2 of the saturated 
magnetization value and is the one proposed by Bos et al 22 . 
The transition is not accompanied by any observable change 
in the broadening of the NMR lines or any clear splitting, 
showing that there is no observable symmetry reduction in 
the crystalline structure from the high-temperature P2i3 space 
group, in agreement with the reported powder XRD dataiS In 
addition, we provide a microscopic description of the prob- 
lem based on a spin Hamiltonian with two nearest-neighbor 
exchange couplings: one antiferromagnetic ( J a f m ~ 68 K) be- 
tween Cui and CU2 ions, and one ferromagnetic ( Jf m ~ —50 
K) between G12 ions. The mean field theory predictions of 
this model are in good qualitative agreement with the behavior 
of the local moments extracted from NMR and magnetization 
data. 

The article is organized as follows. In Sec. [II] we sum- 
marize the synthesis process and give some technical details 
of our experimental procedure. In Sec. IIII Al we present our 
NMR spectra and give a first discussion of some central find- 
ings. The theoretical framework for the explanation of the 
NMR spectra is given in Sec. IIII B I based on symmetry con- 
siderations. The comparison to the NMR data is then given 
in Sec. IIII CI where we also extract the relevant transferred 
hyperfine field parameters as well as the T-dependence of the 
local moments of both types of Cu 2+ ions. A comparison to 
the mean field theory predictions is also made here. Our NMR 
results for the nuclear spin-lattice 1/Ti and spin-spin 1 /T2 re- 
laxation rates are presented in Sec. [IV] We conclude in Sec. 
[V]with a brief discussion of our results. 
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FIG. 2: 77 Se NMR line shape measurements at 220 K in an exter- 
nal field of 14.09 T applied parallel to the [111], [110] and [100] 
crystallographic directions. 



TABLE I: The experimentally observed paramagnetic shifts Av exp 
(in kHz) at T — 220 K measured with the external magnetic field 
Ho along the [111], [110] and [100] crystallographic directions. 



Aexp 
"11101 



A. 



cxp 



Pi 
P2 
P 3 
P4 



-121.8±0.1 
-49.4±0.2 
89.9±0.1 
277.1±0.2 



-142.9 ±0.1 
-53.9 ±0.1 
60.9 ±0.1 
206.8 ±0.1 



-96.6 ±0.01 
137.5 ±0.04 



II. SYNTHESIS AND EXPERIMENTAL DETAILS 

Single crystals of Cu20Se03 were grown by the standard 
chemical vapour phase method. Mixtures of high purity CuO 
(Alfa-Aesar, 99.995%) and Se0 2 (Alfa-Aesar, 99.999%) pow- 
der in molar ratio 2: 1 were sealed in quartz tubes with elec- 
tronic grade HQ as the transport gas for the crystal growth. 
The ampoules were then placed horizontally into a tubu- 
lar two-zone furnaces and heated very slowly by 50°C/h to 
600°C. The optimum temperatures at the source and deposi- 
tion zones for the growth of single crystals have been 610°C 
and 500°C, respectively. After six weeks, many dark green, 
almost black Cu20Se03 crystals with a maximum size of 
8x6x3 mm were obtained. 

The orientation of the crystal axes with respect to the 
crystal faces were determined by Laue X-ray back-scattering 
measurements and by single crystal X-ray diffraction. The 
diffraction data are in agreement with previously published 
datai 22 ' 25,26 The crystal structure of Cu20Se03 which be- 
longs to the cubic space group P2i3 is presented in Fig. Q] 
The structure is a three-dimensional array of distorted corner- 
sharing copper tetrahedra. The unit cell consists of 16 Cu 2+ 
ions which belong to two crystallographically different groups 



denoted here by Cui and Cu?, in 4a and 12b sites respec- 
tively. The oxygen atoms form two different types of distorted 
CuOs polyhedra, i.e., trigonal bipyramidal polyhedra for Cui 
sites and square pyramidal polyhedra for CU2 sites. The CuOs 
polyhedra are connected by sharing edges and corners. Sim- 
ilarly, the unit cell contains two crystallographically inequiv- 
alent groups of Se 4+ ions, Sei and Sen, each one with multi- 
plicity 4. The two types of SeC<3 (lone pair) trigonal pyramids 
share corners with the CUO5 polyhedra. 

Pulsed NMR experiments were carried out on 77 Se nuclei 
with nuclear spin I = 1/2, gyromagnetic ratio ^n/2tt = 
8.1179 MHz/T and natural abundance 7.58%. The NMR ex- 
periments were performed on a 125 mg single crystal at an 
Oxford 14.09 T magnet equipped with a variable temperature 
insert utilizing a home built spectrometer. The echo signal 
was produced by a standard Hahn echo pulse sequence with a 
typical tt/2 pulse length of 6 fj,s. The separation between the 
two echo-generating pulses was 20—60 /is depending on the 
experimental conditions. The ' 7 Se NMR spectra were mea- 
sured by Fourier transform (FT) of the half spin-echo signal 
whenever the whole line could be irradiated with one radiofre- 
quency pulse. For the broad lines, the spectrum was obtained 
by plotting the area of the echo as a function of the irradiation 



frequency (frequency sweep). The nuclear spin-lattice T\ and 
spin-spin T 2 relaxation times were measured using the stan- 
dard spin echo pulse sequence combined with the saturation 
recovery method for T\ measurements. For the analysis of the 
NMR data we have measured the magnetization of a 28.86 mg 
single crystal of Cu20SeC>3 at 14 T using a Quantum Design 
PPMS (Physical Properties Measurement System) apparatus 
located in the Trinity College in Dublin, Ireland. 



III. 77 SE NUCLEAR MAGNETIC RESONANCE SPECTRA 
A. General results 

The nuclear 77 Se spins provide a powerful local probe of 
the behavior of the electronic Cu 2+ moments through the 
transferred hyperfine and the magnetic dipolar mechanisms. 
Here we have performed detailed 77 Se NMR lineshape mea- 
surements at 14.09 T with the magnetic field applied parallel 
to the crystallographic directions [111], [110] and [100]. Fig- 
ure |2] shows the 77 Se NMR spectrum at T = 220 K. There 
appear four distinct resonance lines when the field is applied 
along [111] and [110], while two lines are observed along the 
[100] direction. The specific values of the shifts of the lines 
A;/ = v — vl from the bare Larmor frequency ul are provided 
in Table U The integrated areas under these lines are found 
with the relative ratios 3 : 1 :3 : 1 f or [ 1 1 1 ], 2:2:2:2 for [ 1 10] , and 
4:4 for [100]. Since the area under a given line is proportional 
to the number of the nuclei that resonate in the corresponding 
frequency window, we conclude that there are four magneti- 
cally inequivalent groups of 77 Se sites in the [111] and [110] 
directions and only two groups for [100]. It is clear that the 
latter correspond to the two crystallographically inequivalent 
groups of selenium sites Sei and Sen which have the same 
multiplicity (four per unit cell). On the other hand, both Sei 
and Sen groups split into two magnetically inequivalent sub- 
groups when the field is along [111] and [110] with multiplic- 
ities 3:1 and 2:2 respectively. Below in Sec. IIIIBI we shall 
be able to identify these specific subgroups of Se sites (and 
even retrieve the values of the corresponding elements of the 
transferred hyperfine tensor) by taking into account the local 
symmetry around each Se site under the conditions that (1) we 
are in the 3up-ldown state and (2) that the crystalline structure 
belongs to P2i3 space group. The second condition is an im- 
portant issue since any symmetry reduction of the crystal from 
P2i3 to either of the two possible crystallographic subgroups 
R3 or P2i2i2ir2 will result in partially or fully splitting of 
the above multiple lines. 

Let us now examine what happens as we cool down the 
system. We have studied the temperature dependence of the 
NMR spectra in the range 20-290 K with the magnetic field 
applied along [111] and [110]. Some representative NMR 
spectra are shown in Fig. [3] We first note that the lines PI 
and P2 get gradually closer to each other below 100 K and 
merge into a single peak at lower temperatures. We also find a 
gradual increase in the line broadenings or Full Width at Half 
Maximum (FWHM). However, as we show in Fig. [4] all the 
NMR line broadenings (shown with symbols in Fig. follow 
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FIG. 3: Se NMR line shape measurements for selected temper- 
atures measured at external field of 14.09 T applied parallel to the 
[111] and [110] crystallographic directions. 



quite closely the corresponding gradual increase of the mag- 
netization as measured by PPMS (solid line). This being the 
typical behavior of inhomogeneous broadening, together with 
the fact that we find no observable splitting of any of the mul- 
tiple lines, we are lead to conclude that there is no clear sign 
of any symmetry reduction of the crystalline structure as the 
system enters the ferrimagnetic state. 

Another important issue is the behavior of the line shifts as 
we cool down the system, since this is governed by the corre- 
sponding ordering behavior of the electronic Cu 2+ moments. 
Figure|5]shows the T-dependence of the relative shifts IS.v/vl 
measured at Ho = 14.09 T for fields applied along the direc- 
tions [111] and [110]. From these data, as we show below in 
Sec. IIIIBI we shall be able to extract the T-dependence of the 
local moments of both Cui and Cu2 ions and thus confirm the 
ferrimagnetic nature of the low-T phase of this compound. 



B. Theoretical analysis 

In order to understand and analyze the above NMR data we 
first note that the most relevant terms in the nuclear spin- 1/2 
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FIG. 4: (Color online) Temperature dependence of the 7 'Se NMR 
full width at half maximum (FWHM) for each of the four spectral 
lines measured at 14.09 T along the [111] and [110] crystallographic 
directions. The FWHM for lines PI and P3 along [111] are multi- 
plied by a factor 0.9 and 0.7 respectively in order to match the other 
curves. The temperature dependence of the magnetic susceptibility 
per copper mol (multiplied by a factor of 9.4) as measured with the 
magnetometer at 14 T is shown with a solid line. 



Hamiltonian for the fc-th Se site are (i) the Zeeman coupling 
with the external field Ho, (ii) the transferred-hyperfine in- 
teractions with the relevant neighboring Cu spins and (iii) the 
dipolar interactions with all the electronic Cu spins, namely 



H k = -j N hI k • (1 + K£) • H + H" 



H 



dip 



where K£, denotes 



the chemical shift tensor, while 



(1) 



and H^, lp are the local effective fields resulting from the 
transferred-hyperfine and the dipolar interactions respectively. 

Transferred-hyperfine contribution. — Let us discuss the 
transferred-hyperfine coupling first. The transferred-hyperfine 
interaction is local, i.e. it involves a 77 Se nuclear spin and 
its six nearest neighbor Cu ions, and is mediated through the 
Cu-O-Se bonding. The local transferred-hyperfine field expe- 
rienced by the fc-th nuclear spin can be written as 



(2) 



where q = 1/(7jv%/Ub), Mj = — ^^(Sj) denote the local 
moments of the six neighbor Cu ions, /j,b is Bohr's magne- 
ton, and g is the spectroscopic factor with g ~ 2.11^ To go 
beyond this general description we need to exploit the local 
symmetry properties of the electronic environment of each Se 
site. As we mentioned in the introduction, each unit cell con- 
tains four crystallographically equivalent sites of Sei ions and 
similarly four equivalent sites of Sen ions. We denote these by 
Se™ and Sejj respectively, where n = 1, . . . , 4. We start our 
analysis with the first type of selenium ions. As can be seen in 
Fig-IU eacn Se™ ion sits on a high symmetry crystal site. Each 
oxygen of the trigonal pyramid SeC<3 is connected to two dif- 
ferent types of copper ions, denoted as Cui and Cu2. The 
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FIG. 5: (Color online) Temperature dependence of the relative NMR 
shifts Av/vl at Ho = 14.09 T and for each of the four spectral 
lines observed along the (a) [111] and (b) [110] crystallographic di- 
rections. 



three Sei-0-Cui-Cu2 bonding groups are equivalent and are 
mapped to one another by a rotation of 120° around a 3-fold 
symmetry axis. This C3 axis passes through the selenium site 
and is vertical to the plane of the three Q12 ions. This plane 
is parallel and slightly above the plane formed by the three 
Cui ions, as well as to the plane formed by the three oxygen 
atoms (see Fig. [6]). Thus the transferred hyperfine interaction 
at any given Se™ site can be described by two hyperfine ten- 
sors, ASS and Ag"^ . The former (latter) represents the sum 
of the three hyperfine tensors between the n-th Sei and the 
three Cui (CU2) ions. Now, one of the two 3-fold rotations 
reads [x, y, z) i-> (z,x,y), and this transforms the elements 
of any second rank tensor as 



(3) 



Since this is a symmetry operation we must have xx — yy — 
zz and xy = xz = yz. These conditions describe an axially 
symmetric tensor with the third principal eigenvector along 
the [111] direction. Thus the hyperfine tensors A^ 1 and 
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FIG. 7: (Color online) The local environment of Sen type ions. As in 
the case of Sei ions, there is a 3-fold symmetry axis through the Se 
site which is vertical to the planes of C112 ions. Shown distances are 
in units of A. 



FIG. 6: (Color online) The local environment of each Sei nuclear 
site. There is a 3-fold symmetry axis through the Se site which is 
vertical to the planes of the Cui and CU2 ions. Note that the plane 
of the CU2 ions lies slightly above the plane of the Cui ions. Shown 
distances are in units of A. 



, Cu 2 



are also axially symmetric and their principal axes sys- 
tem can be assigned easily. One is the 3-fold axis (e,3) while 
for the remaining two we may take any pair of mutually per- 
pendicular axes (denoted by ei and e2) on the plane formed 
by the three CU2 (or Cui) neighboring ions. Thus we may 
write 
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where we emphasize that each tensor is written in its own prin- 
cipal axes frame. In what follows the above diagonal elements 
are treated as fitting parameters. 

We now turn to the local symmetry of the SeJ} ions. These 
are bonded through oxygen sites to six Cu 2 ions (Fig. [7} which 
belong to two parallel planes. As in the case of the Se" sites, 
a 3-fold rotation axis is again present and this passes through 
the selenium site of the trigonal pyramid SeC>3 and is vertical 
to the set of planes defined by the Cu 2 ions. The hyperfine 
interactions at each Se 2 site can thus be described with one 
hyperfine tensor, Ag"„ 2 , which stands for the sum of the hy- 
perfine tensors of all six neighboring Cu 2 ions. Due to the 
3-fold axis this tensor is axially symmetric and, as explained 
for the Sei case, the principal axes are the 3-fold axis and the 
two mutually perpendicular axes (their plane is parallel to the 
two Cu planes). Written in its own principal axis system, the 
hyperfine tensor A^?", 2 can thus be written as 



X- 

6' a 




(5) 



Summarizing, the hyperfine tensors for each one of the four 
Se™ sites (respectively SeJJ sites) of the unit cell have the form 
of Eqs. 2] (resp. [5]) as long as they are written in their own 
local principal axis coordinate frame. In Table HI1 we provide 



the coordinates of the unit vector e3 in the (x,y,z) frame for 
each of the eight selenium sites in the unit cell. The coordi- 
nate frames of the hyperfine tensors for the four Se" and the 
four SejJ sites are related with proper rotation of their principal 
axes [ei,e 2 ,e 3 ]. 

Dipolar contribution. — In order to provide a more ac- 
curate quantitative description of the measured NMR spec- 
tra we must also consider the dipolar coupling. In contrast 
to the transferred-hyperfine interactions, the dipolar coupling 
is long-ranged and since we are dealing with a ferrimag- 
net we must also consider the effect of the demagnetization 
field. To this end we follow the Lorentz method- 8 which 
consists of splitting the summation over the dipolar contribu- 
tions into two parts. In the first we perform a discrete sum- 
mation over all individual dipolar contributions from all Cu 
ions enclosed in a sphere of radius R much larger than the 
lattice spacing and much smaller than the size of the sam- 
ple. The summation outside this Lorentz sphere gives rise 
to the demagnetization field H dem which can be evaluated 
using a continuous integration. This gives the well known 
expression^ jj dem = 47r(i — N^) • M, where the demag- 
netizing tensor is specific to the shape of the sample and 
M = (4Mi + 12M 2 ) /a 3 = 4M/a 3 is the total moment di- 
vided by the volume of the sample. The demagnetizing tensor 
is not known but we shall be able to adjust it so that we obtain 
a reasonable T-dependence of the extracted moments at high 
temperatures (cf. below). 

Following the above, the total dipolar field at the fc-th Se 
nuclear site can be written as 



9 £ D i- M ; 
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dem 
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where as usually the cartesian components (a, /3 = x,y, z) of 
the dipolar tensors are given by 



9(Di) =~br 



Ak 
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where tjh stands for the displacement vector from the fc-th 
nuclear site to the j-th electronic spin. Summing over the two 
different types of Cu 2+ ions in Eq. © separately we get 
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TABLE II: The second column gives the unit vector e;i in the fixed 
(x,y,z) frame for each of the eight Se sites per unit cell (first column). 
The following columns give the theory predictions for the sum of 
the dipolar (excluding the demagnetization field) and the transferred 
hyperfine contributions to the relative NMR shifts, when the field is 
along the [111], [110], and [100] crystallographic directions. The 
dimensionless parameters 5, f3, 5, $ are defined in Eqs. <| 1 01>-d 1 31> . 
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Using the known value a 22,25 ' 26 of the positions of the Cu 2+ 
moments one may evaluate the summations over the Lorentz 
sphere with very good accuracy (we find that the sums con- 
verge already for R ~ 5a). Now, the symmetry of the crystal 
necessitates that the resulting dipolar tensors have the same 
principal axes as the corresponding transferred-hyperfine ten- 
sors. Indeed written in the corresponding local coordinate 
frame we find 
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-3.37, 



-0.58, j? ap - 1.17, <5 d 



dip 



where, in units of 10 22 cm 
a dlp - 1.18, /3 dlp = -2.36, «5 dlp 
-1.3, and i? dip = 2.61. 

Predictions for the relative shifts. — Including both the 
transferred-hyperfine and the dipolar contributions it is 
straightforward to deduce the relative shifts for each Se site 
and for each direction of the applied field considered in our 
experiments. The resulting expressions - without the contri- 
bution for the demagnetization field and the chemical shift - 
are provided in Table ITT1 in terms of the dimensionless param- 
eters 

a = (an + a d i p )xi + K r + "di P )X2 (10) 

P = (Ar + /?dip)Xl + (A'r + flfip)X2 (ID 

S = S dlpX i + {5' v + 6' Aip ) X 2 (12) 
= ^dipXi + « + K P )X2 , (B) 

where \i = Mi /Ho and xi = M2/H0 denote the local sus- 
ceptibilities per Cui and CU2 respectively. 



TABLE III: Extracted estimates of the components of the transferred 
hyperfine tensors Ag", 1 , , A^S and Ag^| defined in Eqs. © and $5$ 

(here in units of 10 22 cm~ 3 ) for two possible values of the constant 
ri defined in the text. 
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C. Comparison with experiment 

Let us now compare the theoretical predictions given in Ta- 
ble to our experimental data. We find exact agreement for 
the number of distinct lines as well as for their relative in- 
tensities. When the magnetic field is applied along [111] our 
model predicts that the group of Sei sites give one spectral 
line for Sef and a separate, three times more intense line from 
Sef, Sef and Sef. A similar result holds for the group of Se 
nuclear spins of the second type. In the [110] direction we 
expect four resonance lines with relative ratio 2:2:2:2, while 
in the [100] direction we expect two lines with relative ratio 
4:4 (see Table HI]). These predictions are in perfect agreement 
with our experimental results shown in Figs. [2] and [3] 

For consistency reasons we would like next to contrast the 
T-dependence of a, j3, 8, and -d obtained from the NMR data 
along the [111] direction with that obtained from the NMR 
data along the [110] direction. The results (taken after cor- 
recting the data for the chemical shift and the demagnetization 
field, cf. below) are shown in Fig. [8] and are almost identi- 
cal. This provides a much stronger confirmation of the inter- 
nal consistency of the above theory. In addition, Fig. [8] tells 
us that we may use either the [111] or the [110] data in order 
to extract, in conjunction with Eqs. ([T0l>-(fT3]>, the transferred 
hyperfine field parameters as well as the T-dependence of the 
local moments Mi and M2. In what follows we have taken the 
data along [111]. 

Assignment of the NMR lines. — The next step is to iden- 
tify the specific subgroup of Se sites associated to each given 
NMR line. To this end, we first note that the local fields of 
Se™ and Se{J given in Table [XT] map to one another when in- 
terchanging (a, P) i-> (5, {}). This prevents a straightforward 
assignment of the NMR lines to specific subgroup of selenium 
sites. For example, we cannot infer if it is PI and P2 or rather 
PI and P4 (cf. Fig.0 that come from the same type of Se sites. 
We overcome this drawback by comparing the fits in different 
directions of the applied field. Indeed, by fitting the experi- 
mental data of Ai/ff^j, (see Table Q} with the theoretical ex- 
pressions given in TablellTlwe can reproduce the experimental 
data of both Az^f^ and A^f^ only under the condition that 
PI and P2 come from one group of selenium sites while P3 
and P4 come from the other. This result is further supported 
by the nuclear spin-lattice and spin-spin relaxation time mea- 
surements presented in the following section. However it is 
still not possible at this point to tell whether the pair P1-P2 
comes from type-I or type-II Se sites. We have extracted the 
temperature dependence of the local moments using both pos- 
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FIG. 8: (Color online) A comparison of the values of the parameters 
given in Eqs. d 1 01>- d 131) extracted from the NMR line shifts along the 
[111] direction with those extracted from the [110] direction. 



sibilities and we have found that the choice which gives the 
most physically reasonable behavior is the one which assigns 
the pair P1-P2 to type-II Se sites. 

Hyperfine parameters and local susceptibilities. — We are 
now ready to extract the local moments from the NMR data. 
We first write the measured susceptibility per Cu and the mea- 
sured relative shifts K = Av/vl as 



XtotCT) = Xs P in(T) + Xdk 

K{T) = K spin (T) + K' 



(14) 
(15) 



Here Xspin = (xi + 3x2)/4 is the total spin susceptibility 
per Cu, and Xdia, Xw are the T-independent diamagnetic and 
van Vleck contributions respectively. Similarly, K spm (T) in- 
cludes the dipolar and the transferred hyperfine contributions 
which were given above in Table [II] while K stands for the 
T-independent chemical shift. 

Figure |9ja) shows the so-called Clogston-Jaccarino plot^i 
of the four measured shifts K in both [111] and [110] direc- 
tions versus the measured susceptibility per mol Cu in the high 
temperature regime. The linear behavior shows that both xi 
and \2 are proportional to Xspin in this regime. We can ex- 
ploit this linear behavior in order to extract the values of K c 
for each line if we make an estimate of Xdia + Xvv F° r the 
latter we take the upper bound of Xdia + Xw — 1 x 10~ 4 
cm 3 /mol-Cu which is obtained from the comparison of the 
measured magnetization with the mean field theory predic- 
tion (cf. below). We extract the following estimates respec- 
tively for the lines PI, P2, P3 and P4: K c ~ 5.60 x 10~ 4 , 
8.06 x 10~ 4 , 14.2 x 10~ 4 , 15.7 x 10~ 4 along [111], while 
K c ~ 4.68 x 10~ 4 , 6.73 x 10~ 4 , 13.8 x 10~ 4 , 15.1 x 10~ 4 
along [110]. 

After subtracting the above values of Xdia + Xw and -Kchem 
from the bare data, we adopt the following procedure. We 
first extract an estimate for 6[ r by using the data for 8 at the 
lowest available temperature (T = 20 K) and the relation 
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FIG. 9: (Color online) The fractional shift AupJul for each one of 
the four NMR peaks as measured along (a) [111] and (b) [110], vs. 
the susceptibility \ measured with the magnetometer with T being 
an implicit parameter. The dotted lines indicates the points where 
T = 170 K, and 290 K. 



Mi(20A") = riM(20K), where r\ is a constant which mea- 
sures the effect of quantum fluctuations at low temperatures. 
Using Eq. (Qj), the relation M=(Mi +3M 2 ) /4, and the PPMS 
data for M, we may extract the whole temperature dependence 
of Mi and M 2 . In turn, Eqs. ([T0]i-(fT2l can provide the remain- 
ing transferred hyperfine parameters in a straightforward way. 
For comparison, we have used two values of r\\ The first, 
T\ = —2, corresponds to the case without quantum fluctua- 
tion effects, while the second value, r\ = —1.2, corresponds 
to the reduction in the spin length found by Neutron diffrac- 
tion data.— As for the demagnetization factor, it turns out that 
for Nd < 0.6 the extracted Mi (T) is not positive at room tem- 
peratures but saturates to a negative value. However such a 
behavior would be quite unphysical: Although there is a large 
antiferromagnetic exchange field on each Cui site since it has 
6 neighboring Cu 2 sites, one still expects that Mi should ulti- 
mately change sign at some temperature T* which is above T c 
but certainly well below room temperature. On this issue, the 
mean-field theory described below predicts that T* ~ 152 K. 
So we adjust Nd ~ 0.7, which in conjunction with the value 
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FIG. 10: (Color online) Comparison between the total magnetization 
measured by PPMS (corrected for the van Vleck and diamagnetic 
contributions) and the local moments extracted from NMR (with two 
different values of the ratio n , cf . text) with the corresponding mean 
field theory predictions. We also demonstrate how the value of the 
demagnetizing factor Nd affects the sign of Mi at high temperatures. 



of n can provide a physically more reasonable behavior at 
high temperatures. 

The resulting estimates of the hyperfine parameters are pro- 
vided in TablelllTl while the extracted T-dependence of the lo- 
cal moments is shown in Fig.[l0] The major result contained 
in Fig. [TO] is that the NMR data are fully consistent with the 
3up-ldown ferrimagnetic state, namely that the Cu2 moments 
are aligned parallel and the Cui moments are aligned antipar- 
allel to the applied field. Hence NMR on single crystals pro- 
vides a more direct local probe of this state, mapping out the 
whole T-dependence of the local moments. 

Comparison to microscopic model: — It is worthwhile to 
contrast the above NMR findings with a microscopic spin 
model for Cu20Se03. Following the Kanamori-Goodenough 
rules (see also discussion in Ref. 1221) we introduce two ex- 
change parameters: one ferromagnetic (Jf m ) between Cu2- 
Cu2 ions and one antiferromagnetic ( J a f m ) between CU1-CU2 
ions. Looking at the structure one finds that each Cui neigh- 
bors six Cu2 ions while each Cu2 neighbors four Cu2 ions 
and two Cui ions. The self-consistent equations of the corre- 
sponding mean field theory are given by 



mi 



m 2 



tanh 



2T 



(16) 



where mi. 2 = Mi. 2/(5/1/3/2), b = gfiBHo/ks, and the ex- 
change constants are in units of Kelvin. The mean field transi- 
tion temperature is T C MF = \ (\J fm \ + ^/ Jf m + 3 J a 2 fn] 



while 



per Cu at T > T c are given by the expressions 

(,9Mb) 2 2T-2\J fm \ -3J afm 



Yi 

\2 
Xspin 



8k B 



(T 



T C )(T + T C - 
2T — J a f m 



\Ji 



fm 



Sk B (T - 
{9^B? 2T 



T c )(T + T c ~\J fm \) 
-l/2|J fm |-3/2J afm 



8k B (T - T C ){T + T C - I Jf n - 



(18) 
(19) 
(20) 



the predictions for the local and the total spin susceptibility 



Using the exact numerical solution of Eqs. ( [Tol l and ( fTTT i we 
have obtained a quite accurate fit (cf. Fig. QJJ solid black line) 
of the measured magnetization data in the range 110-300 K, 
using g ~ 2.11^ J afm ~ 68 K and J fm ~ -50 K. This fit puts 
an upper bound on the diamagnetic and van Vleck contribu- 
tions to the susceptibility \da + Xw — 1 x 10~ 4 cm 3 /mol-Cu 
which is of the right order of magnitude. Figure [10] shows 
also the solution for the temperature dependence of the local 
magnetizations Mi and M2 (solid red and blue lines) which 
are to be contrasted with the behavior extracted from NMR. 
The agreement is quite satisfactory. Of particular interest is 
the behavior of Mi at high temperatures. As we discussed 
above, the Cui moments remain antiparallel to the field even 
above T c due to the large negative exchange field that is ex- 
erted from the 6 neighboring CU2 moments. The mean field 
theory prediction for the temperature T* at which Mi even- 
tually changes sign is given by T* = | Jaf m + \Jf m \ — 152 
K. 

So the mean field theory provides a semi-quantitative agree- 
ment with the measured magnetization data and captures the 
essential local physics of the problem, being in agreement 
with the picture obtained from NMR. On the other hand, the 
mean-field theory does not provide a good description at tem- 
peratures close to T c since it neglects long-range correlation 
effects. In addition it can overestimate the value of T c , and 
does not capture quantum fluctuation effects which give rise 
to an overall reduction in the length of the local moments at 
low enough temperatures. 



IV. NUCLEAR SPIN-LATTICE AND SPIN-SPIN 
RELAXATION RATES 

The 77 Se spin-lattice and spin-spin relaxation times, T\ and 
T-2 respectively, were measured for the four different peaks of 
the spectrum when the magnetic field is parallel to the [111] 
axis. We first note that the lines PI and P2 have very similar 
T\ and T-2 values over the whole range of T's investigated here 
and the same holds true for the lines P3 and P4. This is yet 
another confirmation that the lines P1-P2 belong to one type 
of Se sites and the lines P3-P4 belong to the second type. 

The recovery of the longitudinal magnetization follows, for 
each peak, a single exponential behavior in time, as expected 
for a I = 1/2 nuclear spin. The values of spin-lattice relax- 
ation rate 1/Ti derived from the fit of the recovery data are 
shown in Fig. [Ola). At T > 170 K, 1/Ti varies weakly with 
temperature for all lines. We may obtain a theoretical fit of 
1 /Xi by considering the Raman scattering of magnons off the 
77 Se nuclei which is expected to be the dominant process (at 
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FIG. 1 1 : (Color online) Temperature dependence of the Se NMR 
(a) spin-lattice T-f and (b) spin-spin T 2 _1 relaxation rates with the 
magnetic field of 14.09 T applied parallel to [111]. The inset shows 
1/Ti vs. 1/T for the line PI and the corresponding fit (solid line) 
with the contribution from Raman scattering of magnons. 



low enough temperatures) that conserves both the total energy 
and the total angular momentum of the nuclear+electron spin 
system. To this end, we shall make the reasonable assumption 
that the relevant lowest spin-wave excitation branch around 
the minimum A (which is set by the external field) has a typi- 
cal parabolic dispersion of the form = A + A(ak) 2 , where 
a is the lattice constant and A is an energy scale of the order of 
the actual exchange couplings in the system which gives the 
curvature around the minimum of the lowest magnon band. 
Following similar arguments with Refs. l30U3lll one obtains 



line) and gives ^4 2 sin 2 #/A 3 ~ 1.6 x 10 7 erg -1 . Taking A/q ~ 
2 x 10 23 cm~ 3 for the line PI, gives X/ks ~ 50 K which is 
of the right order of magnitude. 

Figure [TTT b) shows the T-dependence of the spin-spin re- 
laxation rate I/T2. These results were obtained by fitting the 
decay of the spin echo signal M(2r), after a tt/2 — t — tt 
pulse sequence, with the proper functional form. We have 
found that the decay of the transverse nuclear magnetization 
follows a single exponential decay in the range 80-290 K 
while it deviates from the single exponential law below 60 
K. The irreversible decay M(2r) of the 77 Se echo signal has 
the following contributions. First, the 77 Se nuclear dipole- 
dipole interaction which can be estimated from the calcula- 
tion of the van Vleck second moment^ 2 . Mf taking into ac- 
co unt the natural abundance of 77 Se.— We have found that 
\J A/f e ~ Se rj 85s" 1 with a slight variation among the dif- 
ferent Se sites, which is much smaller than the experimen- 
tal values reported in Fig.QT| The second contribution to the 
spin echo decay is the Redfield term^ which is of the form 
exp(— 2t/T\). From the 1/Ti values reported in Fig. [TTTa) 
it is obvious that the contribution of the Redfield term in the 
spin echo decay is also very small here. 

The third contribution to the spin echo decay is due to fluc- 
tuating dipolar fields between unlike nuclear spins with the 
most dominant being the one among Se and Cu nuclei. De- 
pending on the time scale (or correlation time r c ) of these 
fluctuating dipolar fields we have the following two limit- 
ing regimes: (i) The fast motion regime, characterized by 
7at (5H z )t c <C 1, where the spin echo decay is a single ex- 
ponential, exp(-2r/T 2 ), with 1/T 2 = ^(SH^t^ and 
7^ (5H%) is the second moment of the selenium-copper dipo- 
lar interaction for which a straightforward lattice sum gives 



7.7 x 10 s . The fast motion approximation is 



M Se-Cu 

in our case applicable for T > 80^. (ii) The quasi-static 
regime, characterized by Jn{5H z )t c 3> 1, where the slow 
fluctuations of the selenium-copper dipolar interaction give 
rise to a non-exponential spin echo decay. In the nearly static 
regime the decay goes as exp[— (2T/T2) 3 ], where I/T2 = 
(7^ (5H 2 Z ) /12r c ) 1 / 3 ^ Identifying r c with the spin-lattice re- 
laxation time of Cu nuclei, which at 30 K is measured to be 
about 1ms,— we find I/T2 ~ 862 s _1 in agreement with our 
low-T data. As we noted above the decay of the transverse 
magnetization is not exponential at low T's. A good fit of the 
data can be obtained by both a square exponential decay (half 
Gaussian) and by a third power exponential. The discrimina- 
tion between them is prevented by the weakness of the NMR 
signal at long times. 



1 

tY 



(A/z 



i 2 sin 2 6> 



2 5 ir 3 h\ 



k B T 
A 



dx 



A/k B T e 



1 



(21) 



where z is the number of Cu 2+ ions closest to the nuclear spin, 
A/z stands for the hyperfine coupling with a single electronic 
spin, and 6 is the angle between the quantization axes of the 
electronic and the nuclear spin. Our optimal fit based on a 
numerical evaluation of the integral of Eq. ( f2TT > with A/ks = 
gUsH (H = 14 T) is shown in the inset of Fig. \TT\ a) (solid 



V. SUMMARY 

We have presented extensive 77 Se NMR measurements in 
single crystals of the magnetoelectric ferrimagnet Cu 2 OSe03. 
The analysis of the data has provided a number of central find- 
ings both for the crystalline as well as the magnetic phase of 
this compound. First, the T-dependence of the two types of 
Cu 2+ moments, extracted from the NMR data, is fully con- 
sistent with a phase transition from the high-T paramagnetic 
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phase to a low-T ferrimagnet whereby 3/4 of the Cu 2+ mo- 
ments are aligned parallel and 1/4 antiparallel to the applied 
field. Below the transition temperature we do not observe any 
clear change in the broadening of the NMR lines or any split- 
ting of the NMR lines, which shows that there is no mea- 
surable symmetry reduction in the crystalline structure from 
its high-T space group P2i3. These results are in agreement 
with previous data from magnetization, high-resolution x-ray, 
and Neutron diffraction measurements on powder samples re- 
ported by Bos et alr^ but also from infrared^, and Raman 24 
studies in single crystals of Cu 2 OSe03. 

We have also developed a microscopic spin model with two 
nearest-neighbor exchange interactions: One antiferromag- 
netic ( J a fm — 68 K) between Cui and Cu 2 ions which forces 
them to be antiparallel to each other, and a second ferromag- 
netic interaction between nearest-neighbor Cu 2 ions which 
aligns then in the same direction, giving rise to the above fer- 
rimagnetic state. We have shown that a mean field solution 
of this model provides a very good description of the physics 
of the problem and is in excellent agreement with measured 
magnetization data in a wide temperature range and, more im- 



portantly, it is also consistent with the local picture extracted 
from NMR. A first-principles study of this system may pro- 
vide a more accurate and refined microscopic spin model.— 

More generally, our NMR study provides a strong local 
confirmation that the ferrimagnetic ordering in Cu20Se03 
does not proceed via a spontaneous lattice distortion and thus 
this material provides a unique example of a metrically cubic 
crystal that allows for piezoelectric as well as linear magneto- 
electric and piezomagnetic coupling. 
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